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Preface 


The  purpose  of  this  study  was  to  find  a  method  of 
measuring  the  concentrations  of  charge  carriers  in  n-type 
semiconductors  within  the  initial  depletion  region.  This 
method  was  to  use  capacitance-voltage  measurements  and  the 
etching  of  layers  from  the  semiconductor  surface.  All  the 
necessary  lab  work  was  carried  out  at  the  Avionics 
Laboratory . 

I  wish  to  express  great  appreciation  to  several  people 
of  the  Avionics  Laboratory  who  helped  me  a  great  deal,  such 
as  Lt  David  Elsaesser,  Capt  Michael  Sopko,  Ben  Carroll, 
Robert  Neidhard,  Michael  Curtis,  Mary  Harshbarger,  and  Sgt 
James  Foster.  I  wish  also  to  express  thanks  to  Dr.  Y.K. 

Yeo ,  my  thesis  chairman. 


Gordon  H.  Gainer,  Jr. 
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Abstract 

Methods  were  developed  for  finding  the  charge  carrier 
concentrations  within  the  initial  depletion  region  of  n-type 
semiconductors.  These  methods  combine  the  capacitance- 
voltage  (C-V)  technique  with  the  etching  of  layers  from  the 
semiconductor  surface.  After  each  etch,  the  depletion 
region  is  made  to  end  at  the  same  location  within  the 
semiconductor.  This  location  is  used  as  a  common  reference 
and  is  arbitrarily  chosen  at  a  distance  deep  below  the 
original  unetched  surface.  For  each  etch  depth,  the  voltage 
which  extends  the  depletion  region  to  the  common  reference 
distance  is  found.  The  voltage  drop  from  an  etch  depth  to 
the  deeper  common  reference  distance  is  the  same  as  it  would 
be  if  the  semiconductor  was  not  etched.  This  voltage  drop 
is  the  sum  of  the  applied  voltage  and  the  built-in 
potential.  The  built-in  potential  depends  on  the  barrier 
potential  at  the  semiconductor  surface  and  the  concentration 
at  the  common  reference  distance.  Assuming  the  barrier 
potential  is  the  same  for  each  etch  surface,  the  built-in 
potential  is  constant  for  each  etch.  Consequently,  the 
applied  bias  voltage  at  each  etch  surface,  can  be  treated  as 
the  potential  at  that  etch  depth.  Then  equations,  such  as 
Poisson's  equation,  can  be  used  to  find  the  carrier 
concentrations  within  the  initial  depletion  region.  These 
methods  were  successfully  applied  to  ideal  computer 


CALCULATING  DEPLETION  REGION  CARRIER  CONCENTRATIONS  WITH 
CAPACITANCE-VOLTAGE  MEASUREMENTS  AND  ETCHING 


1 •  Introduction 

Ion  implantation  of  semiconductors  is  widely  used  in 
fabricating  electronic  and  opto-electronic  devices,  such  as 
FETs  (field  effect  transistors),  ICs  (integrated  circuits), 
and  laser  diodes.  Si-implanted  GaAs  is  widely  used  because 
of  its  high  electron  mobility  and  high  electrical  activation 
efficiency.  Ion  implantation  gives  higher  precision  and 
control  of  dopant  profiles  in  semiconductors  than  epitaxial 
growth  or  thermal  diffusion.  Theoretically,  a  continuous 
variation  of  the  fluence  and  energy  of  the  incident  ions  can 
produce  any  desired  dopant  distribution  (1:9).  In  addition, 
ion  implantation  can  produce  good  uniformity  of  implants 
over  large  regions,  as  well  as  reproducibility  from  wafer  to 
wafer  (2:122).  The  electronic  and  opto-electronic  device 
performances  depend  strongly  upon  the  charge  carrier  depth 
profiles.  For  example,  active  layers  in  microwave  FETs  must 
be  precisely  controlled  to  depths  of  about  1  micron  (3:14). 
Therefore,  it  is  very  important  to  be  able  to  accurately 
measure  the  carrier  profiles.  These  carrier  profiles  are 
usually  obtained  by  the  capacitance-voltage  (C-V)  method  or 
by  Hall-effect  measurements.  However,  both  methods  have 
different  advantages  and  limitations. 


The  differential  Hall  method  can  be  used  for  both  n- 
and  p-type  materials,  but  it  is  destructive.  In  this 
method,  a  Hall-effect  measurement  is  made  of  the  sheet 
carrier  concentration,  the  electrically  active  carriers 
below  the  semiconductor  surface  per  unit  surface  area.  Then 
a  thin  layer  is  removed  by  chemical  etching.  Anther 
Hall-effect  measurement  of  the  sheet  carrier  concentration 
can  be  used  to  find  the  carrier  concentration  in  the  removed 
layer.  A  problem  not  usually  taken  into  account  in 
differential  Hall  profiling  is  that  a  surface  depletion 
region,  produced  by  the  filling  of  acceptor-type  surface 
electronic  states,  causes  the  measured  profile  to  be  shifted 
from  the  true  profile  and  contracted.  However,  a  surface 
depletion  correction  can  be  made  to  the  measured  apparent 
Hall  profiles  (4:5070-5075). 

The  C-V  method  uses  a  reversed  biased  Shottky  barrier 
on  the  semiconductor  to  create  a  variable  depletion  width, 
which  can  be  found  from  its  capacitance.  The  density  of 
carriers  at  the  edge  of  the  depletion  width  can  be  found 
from  the  differential  change  of  capacitance  with  bias 
voltage.  The  C-V  method  is  nondestructive.  However,  the 
C-V  method  does  not  give  a  profile  within  the  initial 
depletion  width  at  zero  bias  voltage,  and  this  depletion 
width  often  is  0.2  *im.  Also,  the  C-V  method  can  be 
practically  applied  only  to  n-type  materials. 


The  objective  of  this  thesis  was  to  develop  a  method  of 


calculating  the  carrier  concentrations  within  the  initial 
depletion  region  of  semiconductors  by  combining  C-V 
measurements  with  etching.  The  method  was  to  be  tested  on 
Si-implanted  GaAs .  The  GaAs  samples  used  were  <100> 
oriented  and  usually  Cr-doped  during  crystal  growth.  The 
samples  had  been  ion  implanted  at  room  temperature  at  an  ion 
energy  of  100  or  200  keV.  The  ion  dose  was  8xl012  or  lxlO13 
cm  .  The  samples  were  annealed  at  800  or  850  °C  for  15 
minutes . 


Sequence  of  Presentation 

Chapter  II  provides  background  information  on 
ion  implantation  and  ion  range  statistics.  Chapter  II  also 
reviews  the  theoretical  basis  of  the  C-V  method  and  the 
cause  of  the  initial  depletion  region.  Chapter  III 
describes  the  Si-implanted  GaAs  sample  preparation  and  C-V 
measurements.  Chapter  IV  discusses  the  theory  of  the 
methods  developed  for  finding  the  depletion  region 
concentrations,  starting  with  Poisson's  equation.  Chapter  V 
describes  the  computer  programs  developed  to  analyze  and 
implement  the  methods.  Chapter  VI  shows  and  discusses  the 
performance  of  the  methods.  Chapter  VII  presents  the 
successful  conclusions  of  this  study  and  the  recommendations 
for  further  study. 
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II .  Background 


Ion  Implantation 

The  electrical  and  optical  properties  of  a 
semiconductor  material  can  be  altered  by  doping  --  adding 
impurities.  The  impurities  can  be  added  before  the 
semiconductor  crystallizes  from  the  liquid  state  or  by 
diffusion  or  ion  implantation  afterwards.  Chromium  is  often 
added  during  crystal  growth  to  electrically  compensate  for 
other  trace  impurities.  Thermal  diffusion  is  used  for 
uniform  doping  over  large  crystal  sections.  In  ion 
implantation,  the  dopants  penetrate  into  the  semiconductor 
as  an  energetic  Ion  beam.  Ion  implantation  has  the 
advantages  of  a  wide  choice  of  dopants,  high  accuracy  in  the 
dopant  concentration  and  area  of  doping,  a  well  defined 
dopant  profile,  slight  lateral  spreading,  good 
reproducibility,  and  high  yields  at  low  cost.  Ion 
implantation  can  also  be  used  at  low  temperatures,  and  it 
can  have  good  uniformity  over  large  areas. 

The  Ion  bombardment  has  the  disadvantage  of  causing 
radiation  damage,  such  as  vacancies,  interstitials,  clusters 
of  point  defects,  dislocations  and  dislocation  loops,  and 
stacking  faults,  within  the  crystal  structure.  Also,  some 
of  the  implanted  ions  do  not  go  into  substitutional  lattice 
sites  to  become  electrically  active.  These  problems  can  be 
mostly  corrected  by  annealing  the  crystal.  At  the  high 


4 


annealing  temperatures  (800  -  1000  °C),  dopant  and  substrate 


atoms  can  diffuse  out  of  the  crystal.  This  is  prevented  by 
encapsulating  the  crystal  surface  before  annealing. 

Lindhard .  Sharf f ,  and  Schiott  ( LSS )  Theory 

In  ion  implantation,  the  dopant  ions  suffer  collisions 
with  the  substrate  nuclei  and  electrons  before  stopping. 

The  LSS  theory  takes  into  account  the  electronic  and  nuclear 
stopping  powers  to  calculate  the  penetration  depth.  The 
penetration  depth  depends  on  the  initial  kinetic  energy  of 
the  incident  ion,  the  atomic  number  and  atomic  mass  of  the 
incident  ion  and  substrate  atoms,  the  temperature,  etc. 
Lindhard,  Sharf f,  and  Schiott  used  the  Thomas-Fermi 
interatomic  potential  to  calculate  the  ion  energy  loss  due 
to  nuclear  collisions  and  the  range  or  depth  that  ions  would 
go  below  the  surface  of  the  substrate.  Assuming  a  Gaussian 
range  distribution,  they  also  calculated  a  mean  square 
fluctuation  in  this  range  (1:36-39).  The  ion  implanted 
dopant  concentration  is  then  given  by 


N(xp)  = 


<P 


0p\f7r 


EXP 


-(x  -R  ) 

P^P 
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2  o. 


(2.1) 


where  N(Xp)  is  the  dopant  concentration  at  a  depth  Xp 
measured  along  the  ion  beam  direction,  <t>  is  the  ion  beam 
fluence,  and  Op  is  the  standard  deviation  in  the  projected 
range  R^.  This  is  illustrated  in  Figure  1.  Gibbons  et  al 
(5)  developed  a  computer  program  to  calculate  the  range  and 
standard  deviation  for  various  implanted  ions  and 


SURFACE 


Figure  1.  LSS  Theoretical  Ion  Depth  Profile 


substrates.  For  example,  Si  ions  with  an  energy  of  100  keV 
have  a  projected  range  in  GaAs  of  0.085  urn  and  a  standard 
deviation  of  0.0442  um. 

The  Gaussian  distribution  is  for  the  random  collisions 
within  an  amorphous  solid.  However,  in  a  crystal,  the  ions 
can  be  channeled  down  a  major  crystal  axis  and  penetrate  to 
much  greater  depths  than  in  an  amorphous  solid.  This 
channeling  is  prevented  by  aligning  the  ion  beam  at  an  angle 
to  the  major  crystal  axis,  as  shown  in  Figure  2  (6:13). 


<100>  Axis 


In  Bern 
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Figure  2.  Ion  Beam  Orientation  to  Prevent  Channelin 


Depletion  Reaion  of  Metal-Semiconductor  Contacts 


When  a  metal  is  brought  into  contact  with  an  n-type 
semiconductor,  the  Fermi  energy  level  Ep  stays  constant 
across  the  metal-semiconductor  boundary  as  shown  in  the 
energy-band  diagram  of  Figure  3.  The  conduction  band 
minimum  Eq  of  the  semiconductor  starts  out  level  at  a 
distance  W  from  the  metal  and  then  rises  to  a  barrier  height 
q#B  at  the  metal  surface,  where  q  is  the  electronic  charge 
and  is  the  barrier  potential.  The  built-in  potential  V 
is  the  difference  between  the  potential  at  the  metal  surface 
and  the  potential  at  W.  This  built-in  potential  transfers 
electrons  to  the  metal  and  ionizes  donor  atoms  from  the 
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metal  surface  to  W,  so  W  Is  called  the  depletion  width. 
Chandra  et  al  (7:646)  showed  that  the  difference  between  the 
conduction  band  minimum  and  the  Fermi  level  is  given  by 


Ec-Ef  =  kT  ln[Nc/N(W)] , 


(2.2) 


where  k  is  the  Boltzmann  constant,  T  is  the  temperature,  Nc 
is  the  density  of  states  at  the  conduction  band  minimum,  and 
N(W)  is  the  carrier  concentration  at  the  depletion  width  W. 
The  carrier  concentration  is 

N=ND-NA-n+p,  (2.3) 

where  is  the  donor  impurity  concentration,  is  the 
acceptor  concentration,  n  is  the  electron  concentration,  and 
p  is  the  hole  concentration.  The  electron  and  hole 
concentrations  are  negligible  where  -  N»  is  high.  The 


built-in  potential  is  then 


Vbi  =  9B  -  (kT/q)ln[Nc/(ND-NA)] . 


(2.4 


According  to  Shottky  (6:4),  the  built-in  potential  Vbi 
should  be  the  difference  between  the  work  functions  of  the 
metal  and  the  semiconductor,  Vbi  «  -  #sc.  However, 

experiments  showed  that  the  rectification  of  metal- 
semiconductor  contacts  can  be  independent  of  the  metal  work 
function.  In  1947,  Bardeen  (9:5)  showed  that  the  surface 
barrier  potentia  1  0g  IS  determined  by  surface  states.  The 
surface  states  produce  the  potential  even  without  a  metal 
contact.  As  explained  by  McKelvey  (10:485),  the  Soviet 
physicist  Tamm  showed  in  1932  that  the  surface  states  are 


caused  by  the  ending  of  a  periodic  potential,  such  as  the 
Kronig-Penney  square  well  potential,  by  a  surface  potential 
barrier,  as  shown  in  Figures  4  and  5.  These  additional 
states  within  the  forbidden  energy  band  of  the  Kronig-Penney 
model  are  the  discrete  allowed  energy  levels  of  wave 
functions  localized  near  the  surface.  Shockley  (11) 
calculated  that  there  is  one  surface  state  for  each  surface 
atom.  Surface  states  can  also  be  created  by  impurity  atoms, 
oxide  layers,  and  surface  imperfections.  Treating  the 
surface  electronic  states  as  unfilled  orbitals  or  dangling 
bonds,  Massies  et  al  (12:64)  experimentally  observed  the 
electron  transition  of  a  Ga  3d  orbital  to  a  dangling  Ga  bond 
on  the  (100)  face  of  GaAs.  The  charging  of  the  surface 
states  causes  the  conduction  band  minimum  and  the  valence 
band  maximum  to  rise  as  they  approach  the  surface,  as  shown 
in  Figure  3.  Henisch  (13:184)  explains  that  the  Fermi  level 
also  rises  to  the  surface,  but  when  the  semiconductor  is 
brought  to  the  metal,  its  Fermi  level  falls  to  equal  the 
Fermi  level  of  the  metal.  This  creates  an  electric  field 
between  the  metal  and  the  semiconductor.  If  the  density  of 
surface  states  is  large  enough  to  take  any  additional 
charges  that  would  be  caused  by  the  electric  field,  without 
much  changing  the  Fermi  level,  then  the  space  charge  in  the 
semiconductor  is  unaffected  by  the  metal  contact,  and  the 
barrier  height  potential  is  independent  of  the  metal  work 


function  (13:184). 


When  a  metal  and  an  n-type  semiconductor  are  brought 
together,  a  space  charge  or  depletion  region  is  formed  as 
shown  in  Figure  3.  When  a  reverse  bias  voltage  is  applied 
to  the  contact ,  more  electrons  leave  the  semiconductor  to 
the  metal,  extending  the  depletion  width  W  deeper  into  the 
semiconductor  as  shown  in  Figure  6.  In  Figure  6,  EFM  is  the 
Fermi  level  in  the  metal,  and  EpS  is  the  Fermi  level  at  the 
depth  W  within  the  semiconductor .  The  space  charge  region 
can  be  treated  as  a  voltage  dependent  parallel  plate 
capacitor,  with  capacitance 


C  =  £A/W, 


(2.5) 


where  C  is  the  capacitance  in  Farads,  £  is  the  permittivity 
of  the  material  in  Farads/cm,  A  is  the  area  of  the  metal 
contact  in  cm2,  and  W  is  the  depletion  width  in  cm.  As  the 
applied  bias  voltage  V  increases,  W  increases,  and  C 
decreases . 

In  1942,  Schottky  had  the  idea  of  using  capacitance 
measurements  to  obtain  dopant  profiles  (14).  Copeland  (15), 
Miller  (16),  and  Nakhamanson  (17)  improved  upon  Schottky' s 
idea  to  obtain  the  following  equations  for  dopant  profiles: 

-C3 

N ( x )  =  - - -  at  x  =  £A/C,  (2.6) 

£  qAz ( dC/dV ) 

where  N(x)  is  the  carrier  concentration  at  depth  x  below  the 
semiconductor  surface.  The  bias  voltage  is  increased  in 
small  increments,  and  the  capacitance  is  measured  for  each 


voltage  to  obtain  dC/dV  at  each  capacitance. 


To  measure  the  capacitance,  the  bias  voltages  are 
actually  superimposed  upon  an  alternating  voltage  at 
frequency  u.  The  measured  capacitance  is 

C*  =  C/(l-w2r2C2) ,  (2.7) 

where  C  is  the  actual  depletion  layer  capacitance,  and  r  is 
a  series  resistance  within  the  semiconductor  (16:1106).  urC 
is  made  small  so  that  C'  is  approximately  C.  This  is  done 
by  reducing  the  metal  contact  area  to  reduce  C,  requiring  r 
to  be  small  in  the  semiconductor ,  and  making  u  as  small  as 
possible  without  ruining  the  signal- to-noise  ratio.  The 
frequency  must  also  be  high  compared  to  the  relaxation  times 
of  traps  in  the  semiconductor  (16:1106).  Traps  hold  on  to 
electrons  much  longer  than  most  dopant  atoms  and  distort  the 
measured  profile.  The  best  frequency  for  GaAs  is  usually 
about  1MHz . 

This  C-V  profiling  method  has  the  advantages  that  1)  it 
measures  only  the  electrically  active  impurity 
concentration,  2)  it  is  completely  automated,  and  3)  it  is 
non-destructive.  The  C-V  method  has  the  disadvantages  that 
1)  it  cannot  be  used  for  concentrations  above  about  1018 
cm'3,  because  voltage  breakdowns  and  current  leaks  occur,  2) 
it  is  practically  restricted  to  n-type  semiconductors,  and 
3)  it  cannot  measure  the  concentrations  within  the  initial 
depletion  region,  which  occurs  at  zero  applied  bias  voltage. 
The  initial  depletion  region  often  ends  0.2</m  or  deeper  into 
the  semiconductor.  Applying  a  forward  bias  voltage  would 
enable  concentration  measurements  closer  to  the 


semiconductor  surface,  except  that  at  even  small  forward 
bias  voltages  a  diffusion  capacitance  becomes  too  great 
(18:94).  This  paper  explains  a  method  of  combining  C-V 
measurements  with  the  etching  of  layers  from  the 
semiconductor  surface  to  measure  the  charge  carrier 
concentrations  within  the  initial  depletion  region. 


III.  Experiment s 


Sample  Preparation 

The  samples  used  in  this  study  were  prepared  and  doped 
before  the  study  began.  The  substrate  materials  used  in  the 
study  were  <100>  oriented  undoped  or  Cr-doped 
semi-insulating  GaAs .  The  samples  were  originally  cut  into 
1/4  inch  squares  from  2  inch  curcular  wafers  of  the  GaAs. 
Prior  to  implantation,  each  sample  was  cleaned  with  basic-H, 
de-ionized  water,  trichloroethylene,  acetone,  and  methanol, 
and  dried  with  blowing  nitrogen  gas.  Then,  the  substrate 
was  free-etched  with  an  H2S04  :  30%H202  :  H20  solution  in  a 
7:1:1  ratio  by  volume  for  3  minutes.  The  samples  were  then 
ion  implanted  with  Si  at  room  temperature  with  an  incident 
energy  of  100  or  200  keV  for  a  dose  of  8xl012  or  lxlO13 
ions/cm-2.  The  ion  beam  was  directed  at  7  degrees  off  the 
< 1 00>  crystal  axis  to  minimize  ion  channeling.  After 
implantation,  the  samples  were  soaked  in  an  HC1  solution  for 
1  minute  to  remove  a  natural  oxide  layer  which  had  grown 
since  implantation.  The  samples  were  then  immediately 

capped  with  a  1000  A  layer  of  Si^N^  using  a  pyrolytic 
deposition  system.  Next,  the  samples  were  annealed  at  800 
or  850  °C  for  15  minutes  in  flowing  hydrogen  gas  to  activate 
the  implanted  ions.  This  produces  high  activation 
efficiencies,  as  shown  by  Kim  (19:27),  who  achieved 
activation  efficiencies  of  at  least  80%  for  6xl012  cm-2 


Si-implants  in  GaAs .  After  annealing,  the  encapsulants  were 
removed  with  a  48%  hydrofluoric  acid  solution,  which 
required  about  5  minutes.  Each  half  inch  square  sample  was 
then  cut  into  4  quarter  inch  square  samples. 

C-V  Measurements  and  Etching 

A  series  of  concentration  profiles  was  measured  for 
each  Si  implanted  GaAs  sample  using  the  C-V  technique.  The 
samples  were  etched  between  the  capacitance  and  voltage 
measurements  to  find  the  carrier  concentrations  within  the 
initial  depletion  region  of  each  sample.  About  ten  etches 
were  performed  on  each  sample  to  acquire  C-V  data  throughout 
the  initial  depletion  region.  The  samples  were  cleaned 
before  each  set  of  capacitance  and  voltage  measurements, 
because  it  was  found  that  the  cleaning  could  very  much 
affect  the  measurements,  especially  if  a  sample  was  not 
cleaned  for  several  days. 

The  C-V  data  was  obtained  using  the  Hewlett-Packard 
Model  4061  Semiconductor /Component  Test  System,  which  was 
controlled  by  a  9845B  H-P  Computer.  A  sample  to  be  profiled 
was  placed  with  the  implanted  face  down  onto  a  suction  chuck 
as  shown  in  Figure  7.  The  suction  chuck  has  three  ducts. 

One  duct  pulls  the  air  out  of  a  circular  groove  on  the  chuck 
creating  a  vacuum  between  the  chuck  and  the  sample. 

Actually,  the  samples  were  too  small  to  cover  the  circular 
groove,  so  a  sheet  of  plastic  had  to  be  placed  over  the  test 
chuck  and  pressed  down  around  the  circular  groove  to  obtain 
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a  vacuum  between  the  test  chuck  and  the  sample.  This  vacuum 
pulls  the  sample  down  to  the  chuck  and  pulls  mercury  up  the 
other  two  ducts.  The  two  mercury  columns  then  make  metal 
contacts  with  the  test  sample.  One  mercury  duct  has  a  much 
larger  contact  area,  and  it  was  connected  to  the  high  side 
of  the  bias  supply.  The  smaller  mercury  duct  was  connected 
to  the  low  side.  Bias  voltages  from  0  to  -5  or  -10  volts 
were  incrementally  applied  to  the  sample.  The  bias  voltages 
were  supplied  by  the  Hewlett-Packard  4140B  DC  Voltage 
Source.  The  4275A  Multi-Frequency  LCR  Meter  supplied  a  0.01 
V,  1  MHz  test  signal  to  measure  the  complex  impedance  and 
determine  the  capacitance  of  an  assumed  equivalent  circuit. 
The  circuit  had  the  capacitance  of  the  space  charge  region, 
or  depletion  region,  in  series  with  the  resistance  of  the 
material  between  the  two  mercury  contacts.  The  capacitance 
of  the  larger  mercury  contact  area  is  much  greater  than  that 
of  the  smaller  mercury  duct  contact  area,  so  the  measured 
capacitance  is  approximately  equal  to  the  actual  depletion 
region  capacitance.  C-V  measurements  were  taken  at  two 
positions  on  each  sample,  and  it  was  found  that  the 
measurements  could  change  somewhat  with  the  location  of  the 
mercury  contact  area.  The  C-V  measurements  were  stored  on 
tape  in  data  files,  which  were  transferred  to  the  AFIT 
computer  system. 

Between  the  C-V  measurements,  etching  was  performed 
with  an  H2S04 : 30*H202 : H20  solution  in  a  1:1:500  ratio  by 
volume,  at  room  temperature  and  with  constant  stirring. 


This  solution  etched  about  80  A  per  minute.  Alongside  the 
sample  to  be  etched,  two  control  samples  were  also  etched. 
The  control  samples  were  half  coated  with  black  wax,  and  the 
black  wax  was  not  removed  until  after  the  last  etching. 

Since  the  control  samples  were  not  etched  under  the  black 
wax,  a  step  height  could  be  measured  on  these  samples  as  a 
total  etched  depth,  and  thus,  the  etch  rate  could  be 
determined.  The  etched  depth  was  measured  using  the 
Sloan-Dektak  IIA  Surface  Profiling  System. 


IV .  Theory 

Three  methods  were  developed  for  determining  the  charge 
carrier  concentration  within  the  initial  depletion  region. 
They  are  the  Charge  Density  Moment  Method,  the  Voltage 
Second  Derivative  Method,  and  the  Voltage  First  Derivative 
Method.  The  Voltage  Second  Derivative  Method  is  the  best 
mathematically ,  but  the  Voltage  First  Derivative  Method  is 
much  less  sensitive  to  experimental  error. 

Charge  Density  Moment  Method 

When  a  metal  is  brought  into  contact  with  an  n-type 
semiconductor  under  a  reverse  bias  applied  voltage  V,  the 
potential  +  varies  with  depth  into  the  semiconductor  as 
shown  in  Figure  8.  The  potential  falls  off  steeply  and  then 
levels  off  at  a  distance  xw.  <P%  is  the  barrier  potential. 

Ec  is  the  conduction  band  minimum  at  xw,  and  EpS  is  the 
Fermi  energy  level  at  x^.  q  is  the  electronic  charge.  V^ 
is  the  built-in  potential.  The  depletion  region  extends 
from  x  =  0  to  x  =  xw .  As  seen  by  equation  (2.4),  the 
built-in  potential  depends  on  the  endpoints  x  =  0  and 
x  =  xw>  The  potential  +  is  governed  by  Poisson's  equation: 


d^+  o  (x )  qN  ( x )  dE 


where  p(x)  is  the  charge  density  at  x,  e  is  the  permittivity 
of  the  semiconductor,  N(x)  is  the  carrier  concentration 


given  by  equation  (2.3)  at  x,  and  E  is  the  electric  field. 


By  integrating  Poisson's  equation  once,  the  electric 
field  is  found  to  be 

E ( x )  =  -fXwN(x)dx  +  Ew,  (4.2) 

{ •  x 

where  Ew  is  the  electric  field  at  xw  (20:209).  This  is  the 
same  result  that  would  be  obtained  using  Gauss'  Law,  since 
the  integral  term  contains  the  total  charge  per  area  from  x 
to  xw.  By  integrating  a  second  time,  the  voltage  drop  from 
the  surface  to  x,,  is  shown  to  be 

V  +  Vbi  =  -  JXwJXwN(x)dxdx  +  XWEW,  (4.3) 

where  V  is  the  applied  reverse  bias  voltage  required  to 
extend  the  depletion  region  to  xw,  and  the  reverse  bias 
voltages  are  considered  positive.  Giacoletto  (20:210) 
pointed  out  that  this  double  integral  can  be  integrated  by 
parts  to  yield 

V  +  Vfe.  =  ^  J*WxN(x)dx  +  xwEw,  (4.4) 

which  relates  the  applied  voltage  to  a  moment  of  the  charge 
density.  If  the  depletion  region  Is  assumed  to  end 
abruptly  as  shown  in  Figure  9,  then  E^  is  zero,  since  there 
are  no  charges  right  of  xw.  Actually,  the  depletion  region 
ends  gradually,  as  shown  in  Figure  10,  due  to  an  electron 
"tail"  at  xw  (18:74-77).  The  electron  concentrations  near 
xw  are  usually  considered  small  and  the  abrupt  approximation 
applied,  so  that  Ew  is  considered  negligible  (21:790).  With 
the  Ew  term  dropped  off,  the  above  equation  is 


q 

V  +  Vhl  =  -  I  xN ( x ) dx .  (4.5) 

1  t  J  0 

From  this  equation,  it  can  be  seen  that  even  if  one  measures 
V,  Vbi,  and  xw,  an  infinite  number  of  N(x)  profiles  could 
give  the  same  result.  There  is  no  unique  N(x)  profile  that 
satisfies  equation  (4.5). 

Giacoletto  (20:215)  also  showed  that  the  voltage  drop 
from  a  position  x0  to  xw  is 

*(x0)  -  +(xw)  =  -  JXw(x-x0)N(x)dx  +  (xw-x0)Ew,  (4.6) 

f  X  Q 

where  xQ  is  between  0  and  xw.  If  xw  is  far  from  both  x=0 
and  x=xQ ,  then  Ew  and  the  shape  of  the  space  charge  region 
(as  shown  in  Figure  10)  should  change  very  little  with 
etching  of  the  semiconductor  between  the  surface  and  xQ . 
Under  these  circumstances,  it  can  be  seen  from  equation 
(4.6),  that  the  potential  drop  from  x0  to  xw  is  the  same 
regardless  of  whether  material  has  been  etched  away  between 
0  and  Xq .  In  that  case, 

V(xQ,xw)  +  Vbi(xQ,xw)  =  +(x0)  -  *(xw),  (4.7) 

where  V(xQ,xw)  and  vbi(x0'xW^  are  resPectivelY  the  applied 
voltage  and  built-in  potential  when  the  depletion  region 
extends  from  an  etch  depth  of  xQ  to  xw. 

If  one  designates  the  ith  etch  thickness  as  t^  and  the 
distance  of  the  ith  etched  surface  from  the  original 
unetched  surface  as  T, ,  then 


Ti  =  2  t j  and  TQ  =  0.  (4.8) 

j=0 

The  potential  difference  between  and  xw  is 

q  fw 

V ( T i , Xw )  +  vbi(Ti'xW^  =  “  I  (x-Ti)N(x)dx  +  (xw-T1 )EW(TA ,xw) , 

Ti 

(4.9) 

where  Ew(T^,xw)  is  the  electric  field  at  x^  when  the 
depletion  region  starts  at  the  etch  depth  T^  ,  and  xw  is 
given  by 


xw  =  eA/C(Ti(xw).  (4.10) 

C(Ti(xw)  is  the  capacitance  measured  when  the  applied  bias 
voltage  is  V(T. ,5^) .  There  are  various  ways  equation  (4.9) 
can  be  used  to  find  the  concentration  N(x)  within  the 
initial  depletion  region.  The  voltage  drop  from  the  etch 
depth  T^ _ 1  to  xw  is 


XW 


v<Ti-l-xw)  +  vbi<Ti-l'xw)  =  ~ 


( x-T^_ j ) N ( x ) dx  + 


Ti-1 


+  (xw  Ti-i  )Ew(Ti-i 'xw)  '  (4.H) 

As  shown  by  equation  (2.4),  the  built-in  potent  lal  vbi 
depends  on  the  etch  surface  barrier  potential  and  the 
carrier  concentration  N(xw)  at  xw-  If  one  assumes  that  the 
surface  barrier  potential  is  the  same  for  each  etch 
depth,  and  if  xw  is  used  as  a  common  reference,  then  the 
built-in  potentials  are  the  same: 


vbi(Ti'xW}  "  vbi <Ti-l 'XW> • 


(4.12) 


■  ~  jf -  ■  -  j« - j» -  w  -  w ->•  -  »  -  '■jr-j*  ’M* ■.* ”Jt "  "J*  'mJF v  V  V.V r ,,’-T7^vjrvj»Til 


As  described  earlier,  the  electric  field  at  xw  changes 
little  with  a  change  in  the  etch  depth,  especially  if  xw  is 
far  from  the  etch  surfaces.  The  electric  field  at  xw  is 
usually  considered  as  near  zero  anyway.  Therefore, 
subtracting  equation  (4.9)  from  equation  (4.11)  yields 


which  can  be  expressed  as 

T  , 
r  ■* 

|  (x-T<_i ) N ( x ) dx  = 


,  ?W 

~ [v ( T  j  _  j , xw )  -  V(Ti#xw)]  -  (Ti-Tj.!)!  N(x’ 


(4.14) 


With  the  Nth  etch  depth  T^  as  the  distance  to  the  end  of  the 
initial  depletion  region  of  the  original  unetched  surface, 
as  shown  in  Figure  11,  then 


( x~tn- l )N(x)dx  = 


(  fw 

~[v(TN_i,xw)  -  V(TN,xw)j  -  ( tn_TN- 1 >  N(x)dx. 
q  T„ 


4.15: 


The  quantities  on  the  right  side  of  equation  (4.15)  are  all 


measurable  or  known.  If  the  etch  thicknesses  are  small,  the 
concentration  N(x)  can  be  treated  as  uniform  or  constant 
between  T™,  and  TN.  Using  tN  for  the  Nth  etch  thickness 


Figure  11.  Etch  Depths  Across  Carrier  Profile 


between  the  etch  depths,  the  integral  on  the  left  side  can 
then  be  evaluated  and  the  equation  solved  for  N  to  give 


Ns—  |^[V(TN_1,XW)-V(TN,XW)]  -  tNJ  N{x)dxj  , 

*§  t  TN  j 


4.16) 


and  this  value  of  N  can  be  assigned  to  the  midpoint  of  the 
two  etch  depths: 


2 


With  N ( x )  known  between  TN-1  and  TN,  the  same  method  can 
then  be  applied  to  find  N(x)  between  Tjj_2  and  Tjj.j.  This 
method  can  be  repeatedly  applied  working  backward  to  the 


original  unetched  surface. 


Since  this  method  treats  the  concentration  as  uniform 
between  each  pair  of  etch  depths,  its  values  of  the 
concentrations  will  be  slightly  different  from  the  true 
values.  If  the  value  calculated  for  the  concentration 
within  one  etch  layer  is  smaller  than  the  true  value,  the 
calculated  value  in  the  next  etch  layer  will  compensate  by 
being  higher  than  its  true  value.  Thus,  the  concentrations 
calculated  by  this  method  oscillate  about  the  true  profile. 
In  this  case,  a  best  fit  curve  should  be  obtained.  Another 
disadvantage  of  this  method  is  that  it  cannot  be  applied 
until  after  etching  to  the  end  of  the  initial  depletion 
region . 

Voltage  Second  Derivative  Method 

As  explained  in  the  previous  section,  the  basic  shape 
of  the  space  charge  region  should  stay  nearly  constant  from 
one  etch  depth  to  the  next,  if  xw  i9  far  from  the  etch 
depths  and  the  etch  thickness  is  small.  Then,  equation 
(4.7)  can  be  used,  and 

V(Ti(xw)  =  +(T1)  -  *(xw)  -  Vbl(Ti(xw).  (4.18) 

Furthermore,  if  xw  is  used  as  a  common  reference  distance 
for  each  etch  depth,  *(x^)  and  vt>i^i,xW^  can  treated  as 
constants.  Then,  the  applied  voltage  is  the  potential 
difference  above  +(xw)  +  Vb^ ,  as  shown  in  Figure  12. 
Consequently,  equation  i4.1)  can  be  used  to  find  the 


concentration : 

N  ( x ) 

-t  d2V 
q  dx2 

*  (x) 

1 

[v<T0.xw> 

i 

! 

1 

1 

J 

VV(Tj,xw) 

Vv(T2,xw) 

Figure  12.  Potential  Versus  Etch  Depth 


Considering  the  reverse  bias  voltages  as  positive,  and 
expressing  the  voltage  second  derivative  in  terms  of  etch 
thicknesses  and  etch  depths,  the  concentration  is 


VITj., ,xu)-V(T 


V(Ti  ,  Xrj  )  -  V  (  T  J  ,  1  ,  Xrj  ) 


(4.20' 

and  this  value  of  N  can  be  assigned  to  the  midpoint  of  the 
two  etch  layers: 


It  can  be  mathematically  shown  that  equation  (4.20)  gives  a 
value  of  the  concentration  which  is  between  the  minimum  and 
maximum  values  of  N(x)  in  the  range  between  and  T^  +  1, 

if  the  measured  voltages  and  thicknesses  are  exact. 

However,  since  equation  (4.20)  is  based  on  a  second 
derivative,  it  is  very  sensitive  to  experimental  error. 

Much  higher  accuracy  can  be  obtained  by  smoothing  the  curve 
of  the  voltage  first  derivative  and  then  taking  the 
derivative  of  this  curve  to  obtain  the  voltage  second 
derivative.  Further  accuracy  can  be  obtained  by  using  small 
etch  thicknesses  and  averaging  N(x)  over  a  range  of  x^'s. 
This  Voltage  Second  Derivative  Method  is  much  simpler  than 
the  Charge  Density  Moment  Method  and  can  be  immediately 
applied  after  the  first  two  etches  (three,  counting  the 
zeroth  etch) . 


Voltage  First  Derivative  Method 

By  Gauss'  Law,  the  electric  field  at  a  depth  within  the 
semiconductor  is  proportional  to  the  charge  per  area  past 
that  depth: 


E(x)  “  —  J  N  ( x )  dx  , 


4.22 


If  one  assumes  that  most  of  the  charge  is  within  the  the 


depletion  width,  then 


q  n 

' x )  =  -  [  N(x)dx, 
(  J 


4.23) 


.  V. 


where  x^  is  the  end  of  the  depletion  region.  The  average 


value  of  N ( x )  between  x  and  xw  is 


Nav(x,xw)  = 


fXwN(x)dx 
J  x  t 


E(x) 


(xw-x) 


q  (xM-x) 


4.24 


As  illustrated  in  Figure  12,  the  derivative  of  the  potential 
is  the  derivative  of  the  applied  voltage  which  extends  the 
depletion  region  from  an  etch  depth  x  to  xw.  The  electric 


field  at  x  is  then 


d*  dV 

E  ( x )  - - - - . 

dx  dx 


(4.25) 


In  terms  of  etch  depths  and  thicknesses,  the  electric  field 


fvtT^,  ,xu)-V(T<  ,xu)l 


Substituting  the  electric  field  of  equation  (4.26)  into 


equation  (4.24)  gives 


T  .  +T . 


[v(T1,xw)-V(Ti+1 ,xw)J 


av  2  '  W  _  g 


(4.27; 


*1+1  xw  - 


T  a  +T  . 


This  average  value  of  N  can  be  assigned  to  the  midpoint  of 
the  range  of  its  average: 


•p  +T  ^ 

X  =  - 1±1  +  xw  /2  =  (Td  +  +  xw)/2  .  (4.28) 


This  is  the  Voltage  First  Derivative  Method.  It  is  a  crude 


approximation,  and  it  approximates  the  concentrations  only 


in  the  seconi  half  of  the  depletion  region,  but  it  is  much 
less  sensitive  to  experimental  error  than  the  Voltage 
Second  Derivative  Method.  Also,  the  value  of  N  throughout 
the  initial  depletion  region  can  be  calculated  more 
precisely  by  using  the  Voltage  First  Derivative  Method  as  a 
first  step  of  the  Voltage  Second  Derivative  Method,  as 
described  below. 

The  value  of  N  between  etch  layers  can  be  determined  by 
using  sheet-carrier  concentrations.  The  sheet-carrier 
concentration  Ns(x0,xw)  is  the  total  charge  carriers  per 
area  from  xQ  to  xw.  The  sheet-carrier  concentration  is 
given  by 

Ns(x0,xw)  =  Nav(x0,xw) (xw-x0)  .  (4.29) 

The  concentration  between  depths  xQ  and  Xj  can  be  determined 
from 


N 


=  x0+xi 


Nav(xO'x; 


a  Ns(x0,xw)  -  N£(Xl,xM) 


:xi"x0) 


(4.30) 


By  using  equations  (4.29)  and  (4.30),  the  value  of  N  between 
etch  layers  can  be  obtained  from 


V .  Computer  Programs 

Computer  programs  were  written  to  implement  the  methods 
developed  in  the  previous  section  and  use  capacitance- 
voltage  data  to  calculate  the  carrier  concentrations  within 
the  initial  depletion  region.  To  better  analyze  these 
methods,  files  of  ideal  capacitance-voltage  data  were 
computer  generated  for  various  dopant  distributions.  C-V 
equations  (2.6)  could  be  applied  in  reverse  to  create  the 
ideal  data,  but  it  was  found  that  the  C-V  equations  can  be 
derived  from  the  differential  form  of  equation  (4.5): 


dV 


q 

—  xN ( x ) dx 

e 


(5.1) 


by  using 


C  =  €  A/x  .  (5.2) 

These  equations  (5.1)  and  (5.2)  were  used  to  generate  the 
ideal  C-V  data.  Since  equations  (5.1)  and  (5.2)  are  the 
same  type  of  equations  used  in  developing  the  methods  of 
finding  the  initial  depletion  region  concentrations,  the 
methods  could  not  be  proven  using  the  ideal  computer 
generated  C-V  data.  However,  it  is  assuring  that  the 
methods  are  based  on  the  same  assumptions  as  the  regular  C-V 
technique,  and  the  ideal  data  could  still  be  used  in 
analyzing  the  methods  and  the  effects  of  experimental  error. 
The  capacitances  after  the  ith  etch  were  easily  calculated 
as 


si 


V  -»* v ' 

sJa.  ■*--  A.  ■ 
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C  (  i  ,  j )  =  £A/(x  -Ta)  .  (5.3) 

where  C(i,j)  is  the  jth  capacitance  after  the  ith  etch  and 
Xj  is  the  jth  distance  from  the  original  unetched  surface. 

Xj  =  Ta  +  jh,  (5.4) 

where  h  is  a  small  distance  interval.  Similar  to  equation 
(5.1),  the  voltages  were  determined  from 

[N(Xj)  +  N(Xj-h)] 

V(i,  j)  -  V  (  i  ,  j-1  )  =  —  jh  -  h  ,  (5.5) 

£  2 

where  V(i,j)  is  the  jth  total  voltage  after  the  ith  etch, 
considering  reverse  bias  voltages  as  negative.  A  constant 
built-in  potential  of  0.8  V  was  assumed  throughout  the 
calculations  for  the  ideal  C-V  data,  so  the  zeroth  voltage 
was  set  at  0.8:  V(i,0)  =  0.8  .  The  voltage  increm  nts  were 

added  to  this  starting  voltage  to  obtain  the  applied  bias 
voltages.  The  voltage  and  capacitance  pairs  with  voltages 
less  than  or  equal  to  0,  were  stored  in  a  data  file  for  the 
ith  etch.  Usually,  11  data  files  were  made  for  10  etches, 

100  A  each. 

To  calculate  the  concentrations  from  the  data  files  by 
the  regular  C-V  technique,  the  finite  difference  form  of 
equation  (2.6)  was  used: 

[v  (  i  ,  j  —  1 )  -  V  (  i  ,  j  )  ] 

N(x)  =  - 

£  qA2  r C ( i , j )  -  C  C i . j-1 )  1 

L  J 

where 
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As  expected,  the  profiles  calculated  by  this  C-V  technique, 
for  the  different  etches,  overlapped  exactly  for  the  ideal 
C-V  data.  For  the  experimental  data,  the  profiles  for  the 
different  etches  overalapped  only  approximate 1 y ,  as  shown  in 
Figures  26  to  37. 

The  computer  programs  used  to  find  the  depletion  region 
concentrations  all  required  a  measure  of  V(T^,xw),  the 
applied  voltage  necessary  to  extend  the  depletion  region 
from  the  etch  depth  to  the  common  reference  distance  xw. 
The  distances  from  the  original  unetched  surface  to  the  edge 
of  the  depletion  region  were  given  by 

Dli'jl  'T>*^77  ' 

Thus,  for  each  V{i,j),  there  was  an  associated  distance 
D(i,j).  V(T^,xw)  could  then  be  found  through  linear  or 
spline  interpolation  between  these  distances. 

For  the  Voltage  Second  Derivative  Method,  the 
concentration  at  each  depth  was  averaged  over  a  range  of 
common  reference  distances  xw.  These  ranges  were  chosen 
where  the  C-V  profiles  overlapped  best  for  the  different 
etches.  However,  the  Voltage  First  Derivative  Method  is 
most  accurate  when  the  common  reference  distance  is  chosen 
as  close  as  possible  to  the  original  unetched  surface,  since 


this  narrows  the  distances  over  which  the  concentrations  are 


averaged.  For  this  reason,  each  point  calculated  by  the 
Voltage  First  Derivative  Method  was  calculated  using  only 
one  common  reference  distance,  which  was  as  close  as 
possible  to  the  original  unetched  surface. 


VI.  Results  and  Discussion 


Results  for  Ideal  C-V  Data 

Figure  14  shows  the  performance  of  the  Charge  Density 
Moment  Method  for  the  ideal  C-V  data  of  the  LSS  Gaussian 
distribution  shown  in  Figure  13.  The  calculated 
concentrations  in  the  initial  depletion  region  oscillated 
more  than  expected  about  the  true  curve.  The  amount  of 
oscillation  varied  with  the  choice  of  the  common  reference 
distance  x^j.  This  was  probably  because  of  inaccuracies  in 
the  voltage  interpolations  caused  by  too  large  of  a  step 
size  in  generating  the  ideal  data.  The  Charge  Density 
Moment  Method  cannot  be  applied  until  after  etching  to  the 
end  of  the  initial  depletion  region.  Also,  this  method 
required  an  integration  to  find  the  area  under  the  profile 
from  the  edge  of  the  depletion  region  to  the  common 
reference  distance.  Slight  differences  in  this  area  could 
very  much  affect  the  results  of  this  method,  and  the 
profiles  obtained  using  the  regular  C-V  technique  did  not 
overlap  well.  For  these  reasons  and  because  of  the 
oscillations  about  the  true  profile,  this  method  was  not 
attempted  on  experimental  data. 

The  Voltage  Second  Derivative  Method  worked  very  well 
on  the  ideal  data,  as  shown  for  the  LSS  profile  of  Figure 
15,  the  linear  profile  of  Figure  16,  and  the  parabolic 
profile  of  Figure  17.  The  original  LSS  profile  is  shown  in 


Figure  12,  and  the  original  parabolic  profile  of  Figure  17 
is  shown  in  Figure  18.  In  these  figures,  the  first  data 
point  calculated  by  the  regular  C-V  technique  after  each 
etch  is  labeled  with  the  number  of  the  etch.  Figure  19 
shows  the  results  of  the  Voltage  Second  Derivative  Method 
when  the  voltages  of  every  other  etch  are  purposely  offset 
by  1  mV.  This  was  done  to  show  how  sensitive  the  method  is 
to  experimental  inaccuracies . 

The  Voltage  First  Derivative  Method  worked  much  better 
than  expected,  as  shown  for  the  LSS ,  linear,  and  parabolic 
distributions  in  Figures  20,  21,  and  22,  respectively.  The 

concentrations  calculated  by  this  method  start  at  the 
midpoint  of  the  initial  depletion  region  and  extend  outside 
the  depletion  region,  almost  overlapping  the  profiles 
calculated  by  the  regular  C-V  technique. 

Results  for  Experimental  C-V  Data 

As  explained  at  the  end  of  Section  IV,  the  Voltage 
Second  Derivative  Method  is  based  on  the  Voltage  First 
Derivative  Method,  so  inaccuracies  in  the  latter  show  up 
much  greater  in  the  former.  The  results  of  these  two 
methods  should  be  viewed  together.  For  comparison  with 
these  methods,  an  LSS  distribution  is  shown  in  Figure  23. 

The  LSS  curve  is  for  Si  Implanted  into  GaAs  at  an  energy  of 

■i  o  _  n 

100  keV  and  a  dose  of  1x10  cm  ,  which  was  usually  used  in 
this  study.  The  results  of  the  two  methods  are  shown  in 
Figures  24  -  33.  The  Voltage  First  Derivative  Method  gave 


concentrations  starting  at  the  midpoint  of  the  initial 
depletion  region  and  then  approximately  overlapping  the  C-V 
profiles  just  as  demonstrated  in  Figure  20  for  ideal  data. 
However,  the  curves  calculated  by  this  method  were  much  more 
jagged  for  the  experimental  data.  This  jaggedness  can  be 
seen  as  fluctuations  about  a  true  profile,  so  a  smooth  curve 
can  be  drawn  through  these  fluctuations.  By  using  the  data 
of  a  smooth  curve  with  equation  (4.31),  which  is  another 
form  of  the  Voltage  Second  Derivative  Method,  good  results 
can  be  obtained  for  the  concentrations  throughout  the 
initial  depletion  region.  In  this  way  the  Voltage  First 
Derivative  Method  can  be  used  as  a  first  step  of  the  Voltage 
Second  Derivative  Method.  When  the  Voltage  Second 
Derivative  Method  is  applied  alone,  the  results  can  be 
erratic,  as  shown  in  Figures  25,  26,  28,  30,  31,  and  33.  A 
different  common  reference  distance  was  used  for  each  point 
calculated  by  the  Voltage  First  Derivative  Method,  but  the 
same  common  reference  distance  must  be  used  for  each  point, 
when  the  first  derivative  method  is  used  as  a  first  step  in 
the  second  derivative  method.  This  was  tried,  but  the 
calculated  data  points  were  too  few  and  the  fluctuations  too 
large  to  confidently  draw  a  smooth  curve  through  the 
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Figure  16.  Performance  of  The  Voltage  Second  Derivative 
Method  for  a  Linear  Profile 
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Figure  19.  Performance  of  the  Voltage  Second  Derivative 
Method  When  the  Voltages  of  Every  Other  Etch 
Curve  are  Offset  by  1  mV 
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Figure  21.  Performance  of  the  Voltage  First  Derivative 
Method  for  a  Linear  Profile 
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Pigure  25.  Performance  of  the  Voltage  Second  Derivative 

Method  for  Sample  A311  (Averaged  Over  the  First 
0.2  um  of  Each  C-V  Profile) 
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Figure  26.  Performance  of  the  Voltage  Second  Derivative 

Method  for  Sample  A311  (Averaged  Over  the  First 
0.3  7/m  of  Each  C-V  Profile) 
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Figure  27.  Performance  of  the  Voltage  First  Derivative 
Method  for  Sample  A312 
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Figure  28.  Performance  of  the  Voltage  Second  Derivative 

Method  for  Sample  A312  (Averaged  Over  the  First 
0.3  um  of  Each  C-V  Profile) 
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Figure  29.  Performance  of  the  Voltage  First  Derivative 
Method  for  Sample  S021 
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Figure  30.  Performance  of  the  Voltage  Second  Derivative 

Method  tor  Sample  S021  (Averaged  Over  the  First 
0.3«m  of  Each  C-V  Profile) 
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Figure  32.  Performance  of  the  Voltage  First  Derivative 
Method  for  Sample  S022 
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Figure  33.  Performance  of  the  Voltage  Second  Derivative 

Method  for  Sample  S022  (Averaged  Over  a  Range  of 
0.03  to  0.13  vm  from  the  Start  of  Each  C-V 
Profile) 


VII.  Conclusions  and  Recommendations 


Conclusions 

Methods  were  developed  for  finding  the  charge  carrier 
concentrations  within  the  initial  depletion  region  of  n-type 
semiconductors.  These  methods  combine  the 
capacitance-voltage  (C-V)  technique  with  the  etching  of 
layers  from  the  semiconductor  surface.  After  each  etch,  the 
depletion  region  is  made  to  end  at  the  same  location  within 
the  semiconductor.  This  location  is  used  as  a  common 
reference  and  is  arbitrarily  chosen  at  a  distance  deep  below 
the  original  unetched  surface.  For  each  etch  depth,  the 
voltage  which  extends  the  depletion  region  to  the  common 
reference  distance  is  found.  The  voltage  drop  from  an  etch 
depth  to  the  deeper  common  reference  distance  is  the  same  as 
it  would  be  if  the  semiconductor  was  not  etched.  This 
voltage  drop  is  the  sum  of  the  applied  voltage  and  the 
built-in  potential.  The  built-in  potential  depends  on  the 
barrier  potential  at  the  semiconductor  surface  and  the 
concentration  at  the  common  reference  distance.  Assuming 
the  barrier  potential  is  the  same  for  each  etch  surface,  the 
built-in  potential  is  constant  for  each  etch.  Consequently, 
the  applied  bias  voltage  at  each  etch  surface,  can  be 
treated  as  the  potential  at  that  etch  depth.  Then 
equations,  such  as  Poisson's  equation,  can  be  used  to  find 
the  carrier  concentrations  within  the  initial  depletion 
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region.  These  methods  were  successfully  applied  to  ideal 
computer  generated  C-V  data  as  well  as  experimental  C-V 
measurements  of  Si-implanted  GaAs . 

Recommendations 

Probably  the  main  source  of  error  in  these  methods  is 
the  fact  that  the  depletion  region  ends  gradually  instead  of 
abruptly.  Wu  et  al.  (22:319-329)  showed  how  the  shape  of 
the  depletion  region  changes  with  applied  bias  voltage  for 
the  regular  C-V  technique,  and  in  a  companion  paper 
Klopfenstein  and  Wu  (23:329-333)  described  computational 
methods  for  calculating  the  C-V  measurements  that  would 
result  from  specified  carrier  distributions.  These 
computational  methods  should  be  applied  to  check  the 
accuracy  of  the  methods  developed  in  this  thesis.  The  fact 
that  the  depletion  region  changes  its  shape  should  be  less 
serious  for  the  methods  of  this  thesis  than  for  the  regular 
C-V  technique,  since  a  common  reference  distance  is  always 
used  and  each  point  in  the  depletion  region  should 
experience  about  the  same  potential  regardless  of  whether 
material  has  been  etched  before  it.  If  the  common  reference 
distance  is  far  from  the  etch  surfaces  and  the  etch 
thicknesses  are  small,  the  basic  shape  of  the  depletion 
region  should  change  little  from  one  etch  to  the  next.  This 
is  similar  to  an  assumption  for  the  regular  C-V  technique 
that  the  built-in  potential  changes  little  with  increments 
in  the  applied  bias  voltage. 
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Another  major  source  of  error  is  oxidation  of  the 
semiconductor  surface.  The  surface  barrier  height  was 
assumed  to  remain  constant  for  all  the  etch  surfaces,  but 
oxidation  changes  the  barrier  height.  A  change  in  the 
barrier  height  of  1  mV  can  very  much  change  the 
concentration  calculated.  If  the  barrier  height  could  be 
accurately  measured  for  each  etch  surface,  greater  accuracy 
could  be  obtained  and  the  methods  could  be  simplified. 

For  greater  accuracy,  a  much  larger  number  of  finer 
etches  should  be  performed  both  within  and  outside  the 
initial  depletion  region.  Greater  accuracy  in  the 
measurements  of  etch  thicknesses,  voltages,  and  capacitances 
would  help.  In  this  study,  capacitances  were  measured  to 
within  5,  4,  and  sometimes  only  3  significant  digits.  The 
voltages  may  have  been  measured  to  within  only  3  significant 
digits.  The  calculations  were  rather  insensitive  to 
inaccuracies  in  etch  thicknesses,  but  the  thicknesses  were 
measured  with  an  accuracy  of  about  10%. 
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